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The underlying mechanism of the metal-to-insulator transition (MIT) in BaVS3 is investigated, 
using dynamical mean-field theory in combination with density functional theory. It is shown that 
correlation effects are responsible for a strong charge redistribution, which lowers the occupancy of 
the broader Ai g band in favor of the narrower E g bands. This resolves several discrepancies between 
band theory and the experimental findings, such as the observed value of the charge- density wave 
ordering vector associated with the MIT, and the presence of local moments in the metallic phase. 
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The structural, electronic and magnetic properties of 
the vanadium sulphide compound BaVS3 raise several 
puzzling questions 0, Qj- At room temperature, this 
material crystallizes in a hexagonal (P&z/mmc) struc- 
ture [3(, in which straight chains of face-sharing VSg 
octahedra are directed along the c-axis. At T s ~240 K 
the crystal structure transforms into an orthorhombic 
(Cmc2i) structure 0, thereby creating an anisotropy 
in the afo-plane (i.e., perpendicular to the chain direc- 
tion) and a zigzag distortion of the VS3 chains in the 
6c-plane. Additionally, the Hall coefficient changes sign 
from negative to positve at T s Q ■ On further cooling the 
system displays a metal-to-insulator transition (MIT) at 
T MIT ~70 K. Remarkably, this transition is second-order 
and is not accompanied by magnetic ordering. Only be- 
low T x ~ 30 K indications for an incommensurate anti- 
ferromagnetic order exist 0. 

Forro et al. have found that the MIT can be driven 
to T=0 by applying pressure [7|. Further recent ex- 
periments 0, I3 have demonstrated that the MIT is 
in fact associated with a structural transition. These 
studies establish that a commensurate structural mod- 
ulation sets in, corresponding to a reduced wave vector 
q=(l, 0,^)o in the orthorhombic cell. Furthermore, the 
X-ray diffuse scattering experiments of Fagot et al. |i| 
reveal a large fluctuation regime with critical wave vec- 
tor g c =0.5c*(here c* is the reciprocal unit vector along 
the c-axis of the orthorhombic system), extending up to 
170 K into the metallic phase. In the same temperature 
range the Hall coefficient is strongly increasing |2]. This 
regime might be interpreted as a precursor of the charge 
density wave (CDW) instability, reminiscent of the large 
fluctuations in a quasi one-dimensional (ID) metal in 
the vicinity of a Peierls transition. It should be kept 
in mind however that the conduction anisotropy within 
the system is not strongly pronounced (cr c /cr a ^3-4) |l0j . 
making ID interpretations questionable. The "metallic" 
phase above T MIT displays several other unusual proper- 
ties 0, E3 ■ The resistivity is rather high (a few mfkm) 
and metallic-like (dp/dT>0) only above ~150 K. It dis- 



plays a weak minimum at this temperature, below which 
it increases upon further cooling. Most interestingly, this 
phase displays local moments, as revealed by the Curie- 
Weiss form of the magnetic susceptibility. The effective 
moment corresponds approximately to one localized spin- 
1/2 per two V sites. Since the formal valence is V 4+ , 
corresponding to one electron in the 3(i-shell, this can be 
interpreted as the effective localization of half of the elec- 
trons. At T MIT , the susceptibility rapidly drops, and the 
electronic entropy is strongly suppressed |ll| . 

In the hexagonal high-T phase the low-lying V(3d) 
levels consist of an A lg state and two degenerate E g 
states. A further splitting of the degenerate states oc- 
curs in the orthorhombic phase. First-principles calcu- 
lations of the electronic structure of BaVS3, based on 
density functional theory (DFT) in the local (spin) den- 
sity approximat ion (L(S)DA), have been performed in 
Refs. |l2l H?3 . Q|. For both phases, the calculations do 
yield a V(3(i)-S(3p) hybridization which is strong enough 
to account for the weak anisotropy of the transport prop- 
erties. No band-gap opening has been reached within 
L(S)DA. Instead, very narrow E g bands right at the 
Fermi level, and a nearly filled dispersive band with 
mainly A\ g character extending along the c* direction 
have been found. This is consistent with a simple model 
proposed early on by Massenet et al. [15|. However, the 
occupancy of the narrow E g bands found within LDA is 
too low to account for the observed local moment in the 
metallic phase. The nature of the CDW instability is also 
left unexplained by the DFT calculations. Indeed, the 
calculated value of the Fermi wave vector of the broad 
A lg band is found to be 2fc^ DA ~0.94c* 13], while the 
observed wave vector of the instability is g c =0.5c* [jj. 
Therefore, the simple picture of a CDW at q c —2kp asso- 
ciated only with the A\ g band is untenable within LDA. 
It is likely that the E g states also participate in the in- 
stability, but the LDA band structure does not provide 
a Fermi-surface nesting compatible with the experimen- 
tal q value 1]. Hence, ab-initio calculations based on 
L(S)DA are not sufficient to explain the electronic struc- 
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FIG. 1: LDA band structure for BaVS3 in the Cmc2\ struc- 
ture. Also shown are the fatbands (see text) for the A\ g 
(blue/full dark), E g \ (red/full grey) and E g 2 (green/dashed 
grey) orbital of the V atoms. 
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FIG. 2: Angular momentum-resolved LDA-DOS of the 
symmetry-adapted V(3d) states. 

ture of BaVS3. By using static DFT+U schemes a band 
gap was obtained |16(. However, this required to enforce 
magnetic order, hence leaving the question of the mech- 
anism of the transition into the paramagnetic insulator 
unanswered |24j . 

In this letter, we present calculations in the framework 
of dynamical mean field theory (DMFT), using the LDA 
electronic structure as a starting point. On the basis of 
this LDA+DMFT treatment we propose correlation ef- 
fects in a multi-orbital context as an explanation for the 
discrepancies between band theory predictions and ex- 
periments. Specifically, we show that interorbital charge 
transfers occur, driven by the on-site Hund's coupling. 
This lowers the occupancy of the Ai g orbital in favor 
of the Eg's, hence shifting kp(Ai g ) towards lower val- 
ues. From a calculation of the local susceptibilities, we 
demonstrate that local moments are formed in the metal- 
lic phase, due to the low quasiparticle coherence scale 
induced by the strong correlations (particularly for the 
narrow E g bands). 

On Fig. the band structure of BaVSa calculated 
within LDA for the T=100 K orthorhombic (Cmc2 1 ) 
crystal structure is displayed. The calculations were 
performed by using norm-conserving pseudopotentials 



and a mixed basis consisting of plane waves and localized 
functions |X7l . The results are consistent with previous 
work [lil \lA IT3 |. The symmetry-adapted V(3<i)-basis 
{4> m } w as obtained by diagonalising the orbital density 
matrix n M M>~ Skfc /kXV'kl-^OC^HV'k)) where tpy, stands 
for the pseudo crystal wave function for wave vector k 
and band b, and M,M' denote the cubic harmonics for 
1—2. Being directed along the chain direction between 
pairs of V atoms, the Ai g orbital has mainly d z 2 char- 
acter. In contrast, the E g states, linear combinations of 
d yz , d x 2_ y 2 and d z 2 (E g \) as well as d xy and d xz {E g %), 
only weakly hybridize with their surrounding. The or- 
bitals of the remaining e g manifold point mainly towards 
the sulphur atoms, which results in a large energy split- 
ting, leading to a smaller(larger) contribution to the oc- 
cupied (unoccupied) states well below(above) Ef- Hence, 
the e g states do not have a major influence on the es- 
sential physics around the MIT. Note that in the or- 
thorhombic structure (as already in the hexagonal one), 
the unit-cell contains two formula units, with the two V- 
sites equivalent by symmetry. The high-symmetry points 
T-C-Y in the Brillouin Zone (BZ) define a triangle in the 
k z =0 plane, whereas Z-E-T is the analogous shifted tri- 
angle in the k z =0.5c* plane. The T-Z line corresponds 
to the propagation along the c-axis. We have used a 
"fatband" representation associated with the {Ai g , E g } 
orbitals on Fig. ^ m which the width is proportional to 
the amount of orbital character of each band at a given 
/c-point. Thus, one can identify that the narrow bands 
at the Fermi level are associated with the E g orbitals. 
Along T-Z starting at around -1 eV, one observes the 
dispersive band with strong A\ g character that crosses 
Ef close to the edge of the BZ. Hence the T-Z portion 
of that band is almost filled, with 2Uf— 0.94c* as men- 
tioned above. The E g2 electron pocket at the T-point is 
absent for the hexagonal high-T structure (see also [l^). 
Its existence might be related to the hole-like transport 
below Ts as revealed from Hall measurements 0. The 
partial density of states (DOS) for each orbital displayed 
on Fig. [2 show a rather broad A\ g DOS while the E g 's 
yield very narrow peaks right at and above Ep- 

Recently, it has become possible to investigate corre- 
lation effects in a realistic setting by combining LDA 
with DMFT [3 H<|. Starting from the LDA hamilto- 
nian H^^(k) expressed in a localized basis set, many- 
body terms are introduced, leading to a self-energy ma- 
trix Ti mm i which is taken to be local (k-independent) 
but fully frequency-dependent. In the present work, we 
use a simplified implementation of the LDA+DMFT ap- 
proach, which is geared at keeping those physical ingredi- 
ents which are important for the physics of BaVS3 close 
to the MIT. Specifically, we work within an effective 3- 
band model, whereby the LDA electronic structure enters 
via the DOS in the relevant energy window around the 
Fermi level. The effective bands are derived from the 
symmetry-adapted {Ai g ,E g } states, thus non-diagonal 
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FIG. 3: (a) Band fillings at ^=(fc B T) _1 =15 eV" 1 (T=774 K) 
for the effective bands within LDA+DMFT. (b) Correspond- 
ing shift of the Fermi level for the A\ a band (note that 
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FIG. 4: LDA+DMFT spectral data for [7=3.5 eV, U / J=A. 
(a) V(3d) low-energy quasiparticle bands along T-Z in LDA 
(dashed lines) and LDA+DMFT (solid lines) for T=332 K. 
(b,c) integrated spectral function p(u>) for a single formula 
unit of BaVS 3 at T=1160 K (b) and T=332 K (c). 



J =0). Filled symbols: U/J=7, open symbols: U/J=A. 



self-energy terms E m ^ m ' are negligible. A physically ad- 
equate 'empirical downfolding' procedure was used in or- 
der to construct the effective 3-band DOS. First, we took 
care of the 3d 1 character. The location and width (~2.2 
eV) of the energy window were chosen such that the total 
DOS accomodates two electrons below and ten electrons 
above the Fermi level (per two V). Since the e g bands are 
hardly relevant and of minor weight close to Ep, they 
were hybridized with the E g bands as suggested from 
the resolved partial DOS on Fig. El The contribution of 
these new E g bands was substracted from the total DOS 
normalized to a single formula unit of BaVS3 within the 
chosen energy window. The resulting difference was iden- 
tified as the new downfolded A\ a band, since the A lg or- 
bital substantially hybridizes with the S(3p) orbitals. Fi- 
nally, the self-energies S m associated with these effective 
bands are calculated from LDA+DMFT where the self- 
consistency condition is expressed as an integral over the 
effective partial DOS D^ DA (e). The on-site interaction 
matrix was parametrised as: U}} m —U, UjJu m ,=U—2J 

and f/JJi^/= U— 3 J, with U the on-site Coulomb repul- 
sion and J the local Hund's rule coupling. The DMFT 
local impurity problem was solved using the quantum 
Monte-Carlo (QMC) Hirsch-Fye algorithm. Up to 128 
slices in imaginary time r and at most 10 6 sweeps were 
used [13. 

Fig-BK displays our results for the occupancies of each 
orbital, as a function of U. In the absence of a precise 
determination of this parameter from either experiments 
(e.g., photocmission) or theory (constrained LDA meth- 
ods tend to underestimate the screening for metals), we 
varied U over a rather large range of values. The ratio 
U j J was chosen to be fixed, and two series were stud- 
ied: U/J=7 and U/J=A, corresponding to a weaker and 
stronger Hund's coupling, respectively. The orbital occu- 
pancies in our effective 3-band model, at the LDA level 
(i.e., for U=0) read: n(A lff )=0.712, n(E gl )=0.207 and 
n(E g 2)=0.081. The main effect apparent on Fig. ^ is 



that moderate correlations tend to bring the occupan- 
cies of each orbital closer to one another, i.e., to decrease 
the population of the "extended" orbital A lg and to in- 
crease the occupancy of each E g orbital. For strong corre- 
lations, values close to n(Ai g )~n{E gl )+n{E g 2)— 0.5 are 
obtained in the DMFT calculation, corresponding to a 
half-filled band. In the absence of correlations, it pays to 
occupy dominantly the broad A\ g band, which provides 
the largest kinetic energy gain, while in the presence of 
correlations this has to be balanced versus the potential 
energy cost. The Hund's coupling clearly favors such an 
interorbital charge redistribution, as also pointed out re- 
cently in the context of ruthenates jjjjj. For larger U, 
band-narrowing effects reduce the kinetic e nerg y of the 
Ai g band, hence also reducing its occupancy [21j in favor 
of E g . As expected, a Mott insulating state is obtained 
when U is increased beyond a critical value, which is 
found to depend strongly on J (the range of i/'s dis- 
played in Fig. |3 corresponds to the metallic regime for 
each series, on which we focus in this paper). 

Our calculations reveal that the depletion of the A\ g 
band is accompanied by a reduction of the correspond- 
ing kp along the T-Z direction. (Note that the Lut- 
tinger theorem [22j does not apply separately for each 
band but only relates the total Fermi surface volume 
to the total occupancy). While a full determination of 
the quasiparticle (QP) band structure in the interact- 
ing system requires a determination of the real-frequency 
self-energy, we can extract the low-energy expansion of 
this quantity from our QMC calculation, in the form: 
Re£ m (w + i0+)~Re£ ro (0) - 1/Z m ) + ■■■, with Z m 
the QP residue associated with each orbital. The poles 
of the Green's function determine the QP dispersion re- 
lation: det[w k - Z[H^ BA + ReE(O) - fj]]=Q, with (j, the 
chemical potential. Focusing first on the A\ g sheet of 
the Fermi surface, within our diagonal formulation the 
location of the Fermi wave vector in the interacting sys- 
tem is determined by: e^ A (k F )=fi - RcXU ls (0). This 
quantity therefore yields the energy shift of the Ai g band 
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FIG. 5: T-dependent local spin susceptibilities for U=3.5 
eV, according to the normalization x <loc)= So dr(S z {Qi)S z {r)), 
where S z denotes the z-component of the spin operator. 

at the Fermi surface crossing, as compared to LDA. It is 
depicted in Fig. [^Jd as a function of U . On Fig. 0^,, we 
display (for C/=3.5eV and U/J=A) the QP bands that 
cross the Fermi level along T-Z in a narrow energy range 
around ef- The QP bands are obtained by perform- 
ing a perturbative expansion of the pole equation above, 
which yields: w bk =£ m C b mk Z m [e^ A +Re^ m (0)-^} with 
C mk=l(V'kl</ , m)| 2 the LDA orbital weight. From these 
two figures, it is clear that kp(Ai g ) is reduced in com- 
parison to the LDA value, in line with the global charge 
transfer from A\ g to E g . This opens new possibilities for 
the CDW instability, in particular for the nesting condi- 
tion. 

The enhanced population of the narrow E g bands, as 
well as the correlation-induced reduction of its bandwidth 
(see Fig. 0^,) provide an explanation for the local mo- 
ments observed in the metallic phase. To support this, 
we have calculated (Fig. 03) the local susceptibility as- 
sociated with each orbital Xm° c ' = J2qR- e iXm(q, w=0)]. 
For both values of U/J, the susceptibility of the A\ g 
band saturates to a Pauli-like value at low temperatures. 
In contrast, ^( loc ) of the E g orbitals strongly increases 
as T is lowered (except for the low-filled E g2 orbital at 
U/ J=A). This is because the coherence temperature be- 
low which quasiparticles form is much lower for the E g 
orbitals than for the A\ g . Accordingly, our calculation 
of the integrated spectral functions (Fig. EJ),c) reveals a 
strong T-dependence of the E g QP peak. Some differ- 
ences between the two series are clear from Fig. GD For 
U/J=7, the system is very close to the Mott transition. 
Thus the A\ g electrons also act as local moments over 
part of the temperature range, while for U / J=4, the 
T-dependence of the total local susceptibility is almost 
entirely due to the E g \ electrons. Which of the two sit- 
uations is closest to the physics of BaVS3 does require 
further investigations, albeit some experimental indica- 
tions point at the second possibility |23| . 

In conclusion, we have shown that charge redistribu- 
tions due to correlations lower the occupancy of the broad 
A\ g band in favor of the narrow E g s 1 in comparison to 



LDA calculations. This explains the presence of local mo- 
ments in the metallic phase and paves the road towards 
a full understanding of the CDW instability of this ma- 
terial. Orbital-selective experimental probes are highly 
desirable in order to put to a test the conclusions of the 
present work. Among several outstanding questions still 
open are the detailed nature of the insulating phase (in 
particular regarding the partial suppression of local mo- 
ments and the eventual magnetic ordering) as well as the 
remarkable suppression of T MIT under pressure 0, ■ 
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Note added. After the completion of this work, we 
became aware of a recent angle-resolved photoemission 
study of BaVS3 by S. Mitrovic et al. supporting our 
prediction of a reduced kp(Ai g ). 
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